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ABSTRACT 


An experimental investigation was conducted to determine 
whether any relationship exists between the friction and the heat 
transfer coefficients in the region of nucleate boiling with forced 
convection. The tests were carried out with distilled water flow- 
ing through a 3/8 inch diameter tube. The water was heated by 
passing an electric current through the tube wall. 

Measurements of the velocity of the water were talen at 
threc positions along the test section by rneans of pitot tubes. The 
pitot tubes “s2ie also used to measure static pressure drops. Frorn 
the experimental data the heat transfer and friction coefficients 
were computed. These coefficients were determined for a heat 
transfer range from 1.9 to 2.0 Btu/sq in-sec and for three pres- 
sure levels, 265, 115 and 65 peta. 
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It was found that for a constant subcooling in the nucleate 
boiling region the friction cocificient increases with increasing 
heat transfer rate, indicating that the bubbles actually do infiuence 
the momentum exchange. It was also found that the experimentally 
determined values of the friction and heat transfer coefficients 
plot about a line the equation of which is a = Cyr This is the 
relation that one would expect on the basis of Reynolds’ analogy for 
turbulent flow either without a laminar boundary layer, or with 
the laminar boundary layer when the Prandtl number is unity. It 
appears therefore that the effect of the bubbles is a hydrodynamic 


one leading to an increase in heat exchange as well as in raomen- 


turn exchange. 
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I. INTRODUCTION 


Boiling heat transfer is defined as heat transfer from a 
suriace to a liquid when the temperature of the liquid adjacent to 
the surface is sufficiently high to cause vapor formation. The 
temperature of the bulk of the liquid may be either equal to or 
below the saturation temperature. The difference between the 
saturation temmperature and thé bulk teznperature is called 'sub- 
cooling". 

Fig. 1 is a graph typical of the results found for boiling 
heat transfer under the conditions of either forced or free con- 
vection, with or without subcooling (1). The graph is plotted to 
show the relation of the heat transfer rate (Btu/sq in-sec) to the 
temperature of the heating surface. In the range from A to B 
heat is transferred by convection without boiling. When the sur- 
face ternperature reaches a certain value somewhat above tne 
boiling alll of the liquid, the heat transier rate increases sharp- 
ly with increasing surface ternperature until a maximum is 
reached at point C. Beyond point C the heat transfer rate de- 
creases with increasing suriace temperature. It then passes 
through a minimum at point D, and again increases with further 
_ increase of surface temperature. 

The region from A to B is usually called the convection 
region; the region from B to C, the nucleate boiling region; from 
C to D, the partial film boiling region; and beyond D, the com- 
plete film boiling region. 


At or near point EB bubbles of vapor begin to form on the 
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surface. When the liquid is subcooled, the bubbles will grow to a 
certain size and then collapse again on the surface (1). The maxi- 
mum size is greatly influenced by the amount of subcooling, the 
larger bubbles occurring at the lower subcooling (1,2). With in- 
crease of surface temperature in the nucleate boiling region the 
population of the bubbles increases. Beyond point C the bubble 
population becomes so great that the bubbles begin to coalesce and 
to form fairly large vapor masses. These masses remain for a 
time on the surface and actually form a partial film between the 
surface and the liquid. Beyond point D a stable vapor film has 
covered the entire surface (1). When attempts are made in the 
complete film boiling region to exceed the heat transfer rate which 
can be obtained at point C, the surface material usually melts. 
For this reason point C is often called the ''burnout” point. 

In recent engineering designs high heat transfer rates 
have been required. A propellant cooled rocket nozzle is an 
exarnple of this type of design. Information on boiling heat trans- 
fer is needed for the successful development of such equipment. 
The kind of information most useful to the engineer would perhaps 
consist of a correlation of the burnout point as a function of param- 
eters characterizing the liquid and the flow field. Attempts at ob- 
taining such a relation have, however, not been successful so far, 
and it is felt that further study of the mechanism of boiling heat 
transfer is required before satisfactory relations of this type can 


be derived. 
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One may wish to divide the nucicate boiling heat transfer 
problem into three parts: nucleation of bubbles, bubble motion, and 
finally the effect of bubbles on heat transfer. ‘This paper is con- 
cerned with the last of these aspects. 

It has been well established that the bubbles act to increase 
the heat transfer rate between the surface and the liquid (I, et al.). 
There are several mechanisms which one could postulate to explain 
this improvement (1). Cne suggestion, made by Dr. H. S. Tsien 
has been that ~ at least for boiling heat transfer with forced convec- 
tion - the bubbles act much in the same way as roughnesses. It is 
then believed that the improvement in heat transfer is the same as 
that which could be obtained without boiling with a pipe with equiva- 
lent roughness elements. If this were the case, one would expect 
the friction coefficient to increase as the heat transfer coefficient 
increases. It was the purpose of the tests reported herein to in- 


vestigate the relationship between these two coefficients. 
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Il STATEMENT OF PROBLEM 


As previously stated, the purpose of the present investigation 
was to determine any relationship which may exist between the fric- 
tion coefficient (C..) and the heat transfer coefficient (C,5) for nucleate 
boiling heat transfer. Relationships between these two coefficients 
have been developed on the basis of Reynolds’ analogy for several 
cases of pipe flow (3,4). Ik the Prandtl number of the fluid in ques- 
tion is unity this relationship reduces to the simple equation: 

Cr “¢ 

“a H » (1) 
and this equation should hold for rough as well as for smooth pipes (4). 
If the Prandtl number differs from unity one may still expect Eq. (1) 
to be applicable, provided the roughnesses are large cormmpared to the 
laminar bowendary layer. The results of the present experiments 
will be compared to Eq. (1). 

The heat transfer from, and the pressure drop along a tube 
with flowins water, under conditions of nucleate boiling heat trans - 
fer, have been previously measured with tubes of 1 / 4 inch diameter (5). 
These tubes were so small that the formation of bubbles may have ob- 
6tructed a large percentage of the internal area of the tube and in- 
creased the effective velocity of the water. Without additional infor- 
mation on the velocity, therefore, these pressure crop data could 
not be used for determiining the friction and heat transfer coefficients. 

In the present investigation similar experiments were carried 
out with larger tubes, and an atte:rnpt was made to measure the velo- 


city at three locations along the tube by using pitot tubes inserted in 
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the test section. The pitot tubes were used to measure static pres- 
sure as well as velocity. The pitot tubes (or similar static probes) 
were considered necessary for the static pressure measurements, 

Gince it was not known whether wall piezometer orifices would give 


proper, readings in the presence of bubbles. From these measure- 


ments it was then intended to compute the friction and heat transfer 


coefficients. 
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Wi. TEST INSTALLATION 


A. Test sections 

The test section in which the simultaneous measurements of 
heat transfer and friction were made {see Fig. 2) consisted of a 43 
inch length of 3/& inch diameter tube, with a wall thickness of 0.010 
inch. The tube was made of 304 stainiess steel. All heating was 
done electrically by conducting current through the tube wall. Only 
17 3/4 inches of the tube near the downstream.end were used for 
heating, and for measurements. The first part of the tube, corres- 
ponding to about 56 diameters, was provided so that an equilibrium 
velocity profile would be established before any measurements 
would be taken. Three pitot tubes of 5.053 inch outside diarmeter 
were constructed, placed through drijled holes spaced at seven 
inch inteeveie along the heating section, accurately positioned at 
the centerline-of the test section and fastened by silver soldering. 
The pitot tube measuring stations are numbered consecutively in 
the direction of flow. They will be referred to as stations No. l, 
2 and 3 respectively, throughout the remainder of this paper. For 
reinfercement a small copper sleeve was placed around each pitot 
tube at the silver soidered joint. The first tube waa placed at ap- 
proximately five diameters downstream of the beginning of the heat- 
ing section. Two 3/4 inch diameter copper sleeves approximately 
four inches long were silver soldered to the test section to serve 
as a place of attdchment for electrodes. Stainless steel male tube 
fittings on both ends of the heating portion of the test section served 


to attach the test section into the test piping circuit. A thermal 





ay 


iy 
expansion joint with a Teflon O-Ring was incorporated between the 
tube fitting and the copper sleeve at the downstream end of the test 
section. The entrance end cf the test section wac reinforced with 

2 copper sleeve silver soldered to the test section and machined to 
make a rounded entrance. This test section was used to measure 
velocities at three positions and differential pressures across these 
three positions for various heat transfer rates. 

A second test section wae sirnilarly constructed with the ex- 
ception that no pitot tubes were installed. Instead, four 0.031 inch 
diameter holes were drilled 90° apart in the test section at each of 
the positions corresponding to the pitot tube nressure pick uo posi- 
tions on the first test Fecha: These side holes at each position 
were connected by 2 copper piczomcter zing to which a copper tube 
fitting was salad eth The second test section was necessary in 
order to determine the pressure drop caused by the presence of the 
pitot tubes themselves. T his part of the pressure drop had to be 
subtracted from the readings of the pitot static ontitest in order to 
obtain the net préssure drop caused by the tube walls and the bub- 
bles alone. 

Fig. 2 is a photograph of the two test sections showing de- 
tails of the expansion joint. Fig. 3 is a photograph of a test section 
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which has been cut open to show the placement of a pitot tube. 


B. Test Circuit 
A diagram of the test circuit is shown in Fig. 4. Two fifty 
galion stainless steel tanks were used in the test circuit and were 


connected to the test section, through flanges, by 3/4 inch stainless 
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stecl tubing. The flanges were electrically insulated by Micarta 
and Teflon. An auxiliary piping circuit, for circulating the water 
while heating, was constructed of similar tubing and connected 
the reservoir tank with a circulating pump. The reservoir tank 
and the tubing which lead from the reservoir tank to the test sec- 
tion were heat insulated with asbestos sheeting. Compressed ni- 
trogen was used to force distilled water from the reservoir tank 
upward through the test section into the receiver tank. The test 
section itself was mounted vertically. Choice of the blow-down'' 
cycle instead of a continuous pumped cycle, was made in order to 
take advantage of the inherent relative freedom from pressure 
fluctuations which is characteristic of the blow-down cycle. The 
test installation was located in a laboratory of the Heat Transfer 
Section of the Jet Propulsion Laboratory of the California Institute 


of Technology. 
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C. Instrumentation 

1. Pressure Transmitters 

Five separate differential pressure measurements were 
made on the test section. They were the static pressure drops from 
station No. 1 to station No. 2, and from station No. 2 to station No. 
3, and the three velocity heads at stations No. 1, No. 2 and No. 3. 
All of these pressure readings were taken either by the pitot tubes 
inserted in the test section or through the piezometer holes in the 
tube wall. The orifices and the tubing in each case offer great re- 
sistance to flow because of their gmall size. In order to obtain 


satisfactory and rapid response, therefore, it was necessary to 





m° 
use pressure transmitters which required a minimum of volume 
change. A strain gage type pressure transmitter manufactured 
by the Statham Laboratories Company was found suitable for this 
service, and this type of instrument was used for all five pres- 
sure readings. All of the pressure transmitters were calibrated 
by the Instrument Section of the Jet Propulsion Laboratory and 
were found to be accurate to less than 1°/o of the reading at 6 
inches of mercury in linearity and in hysteresis. The transmit- 
ters were connected to the pitot tubes by 1/8 inch copper tubing, 
and the tubing was electrically insulated from the test section by 
Teflon gaskets between the tube fittings. 

For the test section with the wall piezometer holes, the 
pressure drop between stations was also determined by means of 
Statham gases. These particular transmitters were accurate to 
better than 1°/o of the reading at 2 inches of mercury in linearity 


and in hysteresis. 
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The signals from the transmitters were recorded at the 
central recording station on Leeds and Northrup Speedomax, 
Type G, self balancing continuous recording potentiometers, the 
calibration accuracies of which were kept better than 1 / 2°/o of 


fuli scale readings. 


2. Temperature Measurements 

¥Yive different temperatures were measured at the test 
section: the bulk temperature of the water entering and leaving 
the test section, and the temperatures of the tube wall at stations 


No. 1, No. 2 and No. 3. All measurements were made with 
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Chromel-Alumel thermocouples. The reference temperature in 
each case was melting ice. For the measurements of the bulk temp- 
erature the hot junctions of two thermocouples were each sealed in 
thin glass and placed in the water upstream of the inlet and down- 
stream of the outlet of the test section. Two stainless steel screen 
discs were placed just before the outlet thermocouple to insure 
thorough mixing of the water. In this case the thermocouple read- 
ing is Scat to the water average bulk temperature. Fig. 5 shows 
a photograph of a glass sealed thermocouple in its packing gland 
unit. 

For measurements of the wall temperature the hot junctions 
ef the thermocouples were clamped against the tube wall. They 
were electrically insulated from the wall by a 0.0007 inch maxi- 
mum thickness of mica sheet, and they were separated from the 
clamp by a layer of natural Sil-O-Cel (see Figs. 2 and 6). 

- The eléctrical signals from all the thermocouples were re- 
corded on a Leeds and Northrup Speedomax, Type G, self balancing 
16 channel recording potentiometer. The calibration accuracy of 
this unit was better than 1/2°/o of full scale reading. 

Fig. 6 is a photograph of the test section asa installed. The 
flange attachments at the outlet of the test section and at the inlet 
to the Ha Alina teacbhltes, are for the purpose of electrically insulating 
the test section from the rest of the piping circuit. The pitot tube 
stems are shown protruding from the test section. The leads of 
the wail temperature, thermocouple hot jomeetone may also be seen. 
The outlet bulk temperature thermocouple is shown at the top of the 


photograph, downstream of the flange. In the lower background is 


he 
the ice bath for the thermocouple cold junctions. The pressure trans- 
mitters are shown to the right center of the photograph. The elec- 
trodes from the welding generator power supply are located at the 


top and bottom of the heating portion of the test section. 


3. Power Supply and Additional Equipment 

In order to obtain any desired temperature of the test water, 
it was necessary to be able to heat the water in the reservoir (see 
Fig. 4) prior to the test. For this purpose two 5 kilowatt General 
Electric copper immersion heaters were installed in the reservoir 
tank. In the bottom of the reservoir tank was placed a Weston Ther- 
mometer, Model 1221, with a range of 50 - 500° F, so that the 
temperature of the water could be observed during the heating cycle. 

To regulate the pressure in the water reservoir during the 
discharge cycle, a Grove Loader with a Grove Dome Controller 
was used. The compressed nitrogen was supplied from the central 
system of the ‘laboratory. 

The large amount of current for the test section itself was 
obtained from two Hobart Simplified Arc Welders connected in paral- 
lel. Each welder had a maximum current capacity of 600 amperes 
and a voltage rating of 40 volts. The voltage drop across the test 
section was measured by a Weston D. C. Voltmeter, Model 622, 
and the current was measured by a Weston D. C. millivoltmeter, 
Model 622, in conjunction with a 2000 ampere, 50 millivolt shunt. 
Both meters were accurate to 1/2°/o of the full scale reading. The 
shunt resistance was accurate to + 1/2°%o. 

The water mass flow rate was obtained from a 1/4 inch 


throat diameter Venturi flowmeter constructed at the Jet Propulsion 
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Laboratory. The differential pressure between the venturi throat 
static pressure and the total pressure was recorded on a Barton 
Hydraulic Rotary Recorder, Msdel 202, with a range of 0 - 300 
inches of water. The gage and Venturi flowmeter were calibrated 
together ag a unit and the calibration is shown in Fig. 7. 

An Ashcroit Laboratory Test Gage, range 0 - 300 psi, 
calibrated to better than 1/2°/o accuracy of fuil scale reading, 
was used to read the pressure in the heating tank during the test 
run. 

The water mass flow rate was hand controlled by a 3/4 


inch Vogt stainless steel metering valve. 
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IV. TEST PROCEDURE 


Tests were made with both test sections for each of the three 
pressure levels, 265, 115 and 65 psia, for a mass flow rate of ap- 
proximately 0.75 pounds per second and for heat transfer rates 
ranging from 0.0 to a maximumof approximately 2.1 Btu/sq in-sec. 

Prior to a run, distilled water in the reservoir was heated 
to the desired temperature by the immersion heaters. This temp- 
erature was selected to be approximately 100° F below the satura- 
tion value for each pressure. Cccasionally the water was circu- 
lated by the circulating pump to insure uniform temperature. When 
the water was at the predetermined temperature, compressed ni- 
trogen was introduced to the reservoir. The nitrogen pressure was 
maintained at the desired value throughout the test run by means of 
the regulator. The throttle valve was opened and adjusted to pass 
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the predetermined flow rate. The welding generators were then 
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turned on and ‘the power input to the test section was varied in steps. 
During each run, the temperature of the water entering and leaving 
the test section as well as the three outside wall temperature mea- 
surements were recorded. Simultaneously the velocity was mea- 
sured at the three stations and the differential pressure between 
stations was taken. Prior ” emptying the reservoir of water, elec- 
tric power to the test section was turned off so as to avoid burnout. 
The water in the receiver tank was pumped back into the reservoir 
in preparation for the next run. A similar testing procedure was 
used when the second test section was installed, except that velocity 


and wall temperature measurements were omitted. 
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V. DESCRIPTION OF MEASUREMENTS AND COMPUTATIONS 


A. Test Section Measurements 

The physical measurements of the test sections which were 
used in the computations are tabulated in Appendix B. The outside 
diameter (d.) and the wall thickness (w,) were measured at several 
places near the entrance and exit of the test sections and averaged. 


See Appendix A for a complete list of symbols. 


B. Flowmeter Calibration 

The Venturi flowmeter and the associated pressure gage 
were calibrated together for different gage pressures (see Fig. 7) 
by flowing water at a temperature of 170°F from the reservoir tank 
through the Venturi flowmeter into a tank placed on a platform 
scale. The time for a given weight discharge was measured. Since 
the Venturi flowmeter was calibrated for only one water temperature, 


the mass flow rate (m) had to be corrected for changes in water 


temperature, and consequently in density (p). 


C. Heat Balance 

A heat balance was made for most of the test’runs by com- 
paring the rise in dulk temperature (St, ) of the water through the 
test section, as computed from power measurements, with the ac- 
tual rise. The computed temperature rise was obtained from the 
equation: 

P 3.413 on 
“p"mc tsa0e (2) 

where P is the electrical power in watts, m is given in lbs/sec and 


C is the heat capacity of the water in btu/ lb-°F, The actual bulk 


a &e 
temperature rise in the test section was determined by taking the 
difference of the water bulk temperature at the inlet and outlet of 
the test section. Figs. 8, 9 and 10 are graphs of the actual and 
computed bulk temperature rise for test section pressure levels 
of 265, 115 and 65 psia, respectively, for various heat transfer 
rates (q/A) in Btu/sq in-sec. Only for the runs at 265 psia pres- 
sure was there any detectable heat loss to the surrounding atiiog? 
phere. The loss was approximately 4-1 /2°/0. Figures based on 
the actual bulk temperature rise were used in computing heat 


transfer rates. 


D. Wail Temperature 

The temperature of the outside wall of the test section was 
taken at stations No. 1, No. 2 and No. 3. Figs. 11, 12 and 13 
show a graph of these temperatures for various heat transfer rates 
at various pressures, The temperature of the inside wall (t_ 4) ad- 
jacent to the Lee in the test section was determined by subtract- 
ine the temperature drop through the wall (At) from the tempera- 
ture of the outside wall (to) The temperature drop through the 


wall was calculated from the equation (6): 


w w 
= me i t | 3600x144 o 
aty = A [dat —“—- uy 
where q/A is given in Btu/sq in-sec, w, and d. are given in inches 
and K_ is the thermal conductivity in ee _ (see Fig. 14). 
” ft-hr-°F 


This equation was derived by assuming uniform heating throughout 


the tube wail and no heat flow to the outside. 
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E. Detection of Incipient Soiling 

A graph was made of the outside wall temperature at each 
of the three stations (see for example Fig. 15) on semi-log paper 
tor various heat transfer rates. At the start of boiling the slope 
of these curves changes abruptly. The curves were used to de- 
termine which of the tests were actually conducted in the nucleate 
boiling region. As can be seen from the figure, boiling started 
at all three stations in the heat transfer range between 1.0 and 
1.5 Etu/sq in-sec for a pressure level of 115 psia. Eoiling was 
also found to start within the same range for the pressure levels 


of 265 and 65 nsia. 


FEF. Velocity Measurements 

Velocity measurements were taxen at the centerline of 
the test section and graphs of the measurements are shown in 
Figs. 16, 17 and 18. The absolute values of the velocities at sta- 
tions No. 2 ine No. 3, are not considered reliable since they 
were influenced by the wake of the pimeheiee pitot tube or tubes. 
The change of the velocity readings with increasing heat trans- 
fer rate, however, was taken as an indication of the change in 
average velocity at the station in question. This method of com- 
puting the average velocity was based on the assumption that the 
ratio of the velocity at the center of the test section to the aver- 
age velocity remains unaffected by the rate of heat transfer as 
well as by the bubbles on the wall. This assumption is believed 
to be adequate in view of the fact that the correction involved is 


small. The variation in velocity between measuring stations was 


taxen as linear. 
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G. Pressure Drop : 

For the first test section the pressure drop between stations 
was determined by measuring the difference in static pressure be- 
tween the three pitot tubes. Graphs were made of the pressure 
drops for the three pressure levels, 265, 115 and 65 psia, for var- 
ious heat transfer rates, and for various initial water bulk temper- 
atures (t, 3) (see Figs. 19 to 24). These pressure drop measure- 
ments, however, are not yet the required ones since they include 
a loss which is due to the presence of the pitot tubes themselves. 
This loss must be subtracted from the data obtained from the first 
test section. For this purpose similar pressure drop measurements 
were made along a second test section constructed without pitot 
tubes, but with wall piezometer holes. These pressure drop mea- 
surements were made for heat transfer rates in the non-boiling 
region only (eee Figs. 19 to 24), as the wail pressures could not 
be considered’weliable in the boiling region. The difference be- 
tween the two uty of pressure measurements was taken as the pres- 
sure loss caused by the pitot static tubes. For any given pressure 
level this loss was assumed to remain constant over the whole 
range of heat transfer rates investigated. The pressure loss 
caused by the pitot tubes was then subtracted from the measure- 
ments obtained with the first test section. Curves of pressure drop 


versus heat transfer rate adjusted in this manner are shown in 


Figs. 25 to 30. 


H. Friction Coefficient 


The friction coefficient is defined by the equation: 





pez 


G ; (4) 


F re) 


lowe _ 

tal 4 
where Tis the fluid shear on the tube wall and Vis the average 
velocity. Writing a force balance for a fluid flowing steadily in a 


tube this coefficient can be related to the pressure gradient as fol- 


lows: 


L pV? wd' __ dp xa” 
> 7 =-~ = — 1728, (5) 


where p is given in lbs /tt?, V is given in ft/sec, d' is the effective 
inside diameter, g is the acceleration of gravity in ft/sec’, se is 
the corrected pressure gradient along the test section in lbs /in?. 


The average velocity V, and the diameter d', can be related to the 


mass flow by the equation: 
2 
‘ 
re pV =m. (6) 


Substituting this expression into Eq. (5), one obtains for Cy. the 


relation: 
1/2 


—— & m 1728 
Cp - a Sarin ° 7) 
Vp f 
To determine the gradient ($2), the data shown in Figs. 25 


to 30 were used to construct a graph of pressure versus distance 
along the test section, for three constant heat transfer rates. For 
this graph the pressure at station No. 1} (x = 0 position) was taken 
equal to the pressure in the reservoir tank, the pressure at station 
No. 2 (x = 7 inches) was obtained by subtracting the amount of the 
pressure drop between stations No. 1 and No. 2 (Ap, 2) from the 


pressure at station No. 1, and the pressure at station No. 3 
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(x = 14 inches) waa obtained dyvitat: subtracting the pressure 
drop between stations No. 2 and No. 3. Figs. 31 to 33 show the 
pressure distributions aiong the test section for heat transfer rates 
of 1.9, 1.5 and 2.0 Btu/sq in-sec for each of three pressure levels 
in the reservoir tank, 265, 115 and 65 psia. To determine the 
pressure gradient ($2), it was assumed that the graphs‘of Figs. 31 


to 33 were segments of a parabola, whose equation is of the form: 
2 
pzat+bxtcx. (S) 


The equation was written for each of the three positions x = 0, 7 
and 14 inches and solved simultaneously to determine the constants. 


The derivative, 
SP =b + 2ex, (9) 


was used to determine se at the desired value of x. The local water 
density was of course used for p. The value of V was obtained by 
dividing the make flow rate by the initial water density and the cross - 
sectional area of the tube, and by adjusting this value for any velocity 
increase indicated by the pitot tube readings (see Section V-F). The 


value of m was obtained from the flowmeter data. 


J. Heat Transfer Coefficient 
The heat transfer coefficient is defined for steady flow of 


liquid through a tube as: 


144 9 
6 ww eC, (11) 


a 
ecv{ Ci -t,) 


where the heat capacity (C) is given in Btu/lb-°OF. 


The heat transfer rate (q/A) was determined from the relation: 
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z = At. ae Btu/sq in-sec, (12) 


where At. is the observed actual bulk temperature rise in the water 
through the test section (see Figs. 11 to 13), and A is the area in 
sq in of the heating surface of the test section, based on the actual 
inside diameter. The values of p and of V are the same as those 
ueed for the equation for Co (Eq. 7). The temperature difference 
(t.-t) was obtained by subtracting from the value of the outside 
wall temperature the change of temperature through the test section 


wall (see En. 3) and the bulk temperature of the water in the test 


section. 
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VI, PRESENTATION OF RESULTS 


A. Velocity Measurements 

Figs. 16 to 18 show the reeulte of the velocity measurements 
at the three stations along the test section. The measurements were 
taken at pressure levels of 265, 115, and 65 paia, and they covered a 
heat transfer range from 0.0 to 2.1 Btu/sq in-sec. The velocity mea- 
- gurements were constant at station No. 1 for all three pressures. 
This was to be expected considering that station No. 1 was only 1.7 
inches downstream of the point where heating started. The subcooling 
was still large, and under that condition the bubbles were small, as 
was shown by Gunther (2). The flow velocity at stations No. 2 and 
No. 3 increased somewhat with increase in heat transfer. This in- 
crease became noticeable at about 1.0 to 1.5 Btu/sq in-sec, which 
was the heat transfer rate at which boiling was determined to have 
begun (see Section V-E). The flow velocity at station No. 3 was af- 
fected more than that at station No. Z. This result may be explained 
by the difference in subcooling at the two stations. It is also inter- 
esting to note that at high heat transfer rates, the velocity at station 
Mo. 3, became rather sensitive to the bulk temperature of the water 
at the inlet to the test section. This again is believed to be due to 
the reduced subcooling which resulted from higher inlet temperatures. 

As mentioned before, the absolute values oi the measured 
velocity at stations No. 2 and 3 are not considered significant because 
of the eos of the preceding pitot tube or tubes. The change in velo- 
city with boiling is, however, believed to be indicative of the increase 
in the average velocity, and the average velocity was corrected ac- 


cordingly, as discussed in Section V-F. 
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B. Pressure Measurements 

Results of the pressure drop measurements are shown in 
Figs. 19 to 24. These figures show the results for the test section 
with pitot tubes installed and at pressure levels of 265, 115 and 65 
psia. The pressure drop was plotted as a function of the heat trans - 
fer rate, and values of the water bulk temperature at the inlet to 
the test section are indicated. Figs. 19 to 24 also show similar 
results for the test section without pitot tubes. For the test sec- 
tion with pitot tubes installed, the range of heat transfer rates in- 
vestigated was from 0.0 to 2.0 Btu/sgq in-sec, and for the test 
section without pitot tubes, 0.0 to 1.2 Btu/sq in-sec. The pres- 
eure drop was observed to increase with an increasing heat trans- 
fer rate in all cases after boiling started. The increase was preater 
between stations No. 2 and No. 3 than between stations No. 1] and 
No. 2. This result can again be related to the reduction in sub- 
cooling along the test section. The effect of the inlet water bulk 
termperature may be explained similarly. In the non-boiling re- 
gion the pressure drop was observed to decrease slightly with in- 
creasing heat transfer rate. This is probably due to viscosity 
changes near the tube wall during heating. 

Figs. 31, 32 and 33 represent the sarne pressure data as 
discussed in the previous paragraph, except that they have been 
adjusted for the pressure loss causcd by the pitot tubes as ex- 
plained in Section V-G. 

Finally using the measurements discussed in the foregoing, 
the friction and heat transfer coefficients were computed as out- 


lined in Sections V-H and V-J. The results are presented in two 
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graphs, a graph of the friction coefficient (C.) versus the heat 
transfer rate (q/A) for various pressure levels (Fig. 34), anda 
graph showing the relationship between one-half the friction co- 
efficient and the heat transfer coefficient (Fig. 35). The graphs 


will be discussed in the following section. 





_ a. 
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Vil. CONCLUSIONS 


The fact that nucleate boiling improves the heat E-waeios 
rate is well known (1,2, et al.). The purpose of the present in- 
vestigation was to determine the effect of nucleate boiling on fluid 
friction for a forced convection system. In particular it was de- 
sired to determine whether any relationship exists between fric- 
tion and heat transfer for water flowing through atube. For 
the comparison, the heat transfer and friction coefficients were 
chosen rather than the overall pressure drop and heat transfer 
rates. In this manner it was possible to exclude any effects caused 
simply by changes in the flow velocity or in the temperature differ- 
ential belwe ee’ the tube wall and the water. 

Fig. 34 is a graph of the friction coefficient versus the heat 
transfer rate. The points shown were taken at constant subcooling, 


and at an average pressure of 115 psia. The point at which boiling 


~~ 


oe 


be gan has been‘ indicated. It can be seen that with increasing heat 
transfer rate the friction coefficient actually increases. This shows 
that the increase in pressure drop (see Figs. 25 to 30) cannot be 
explained solely by an increase in the effective velocity which might 
have been cauged by an obstruction of the flow passage by the bub- 
bles. The rise in the friction coefficient indicates that the bubbles 
actually influence the momentum exchange. The simultaneous ef- 
fect of the bubbles on the heat transfer can be examined by consid- 
ering a graph of one-half the friction coefficient (-5=) versus the 
heat thanetes coefficient (C,,). Such a graph is shown in Fig. 35. 


The results shown were taken at three different pressure levels, 
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for a heat transfer range from 1.0 to 2.0 Btu/sq in-sec and ior 
various degrees of subcooling. It can be seen that the experimen- 
tal points define as an average curve, a line given approximately 


by the equation: 


7 Ces. (1) 


This, however, is the relationship which is obtained from Rey- 
nolds analogy either for the case where a laminar boundary layer 
at the heating surface is non-existent or for the case where the 
Prandtl number for the fluid is equal to 1.0. The Prandtl number 
of water varied between 1.0 and 2.0 for the range of bulk temper- 
atures investigated. As the Prandtl] number of the fluid was rela- 
tively close to one, the question of whether or not the laminar 
boundary layer existed will have to be left undecided. In view of 
the results shown in Fig. 35, it can, however, be stated that the 
bubbles affectthe heat transfer and the friction in the ''same way", 
within the meaning of Reynolds analogy. It appears, therefore, 
that the effect of the bubbles is a purely hydrodynamic one, acting 
on the turbulent exchange coefficients both for friction and for heat 
transfer. Thus the results of the present experiments seem to be 
in agreement with Tsien's hypothesis which was mentioned in Sec- 
tion I. 

In conclusion, it should be pointed out that the experiments 
coverdonly a rather limited range of Reynolds numbers (140, 000 


to 380, 000) and that the flow velocities were in all cases higher 
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than the average bubble growth rates. At lower velocities and 
different Reynolds numbers, the nucleate boiling may have a dif-~ 
ferent effect. This should certainly be true for boiling under 
free convection conditions, because in that case the bubble growth 


velocities must have a strong influence on the flow field. 


: ~ 
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APPENDIX A 


Nomenclature 


ifeating area of test section based on length 
and inside diameter of heating section 


Heat capacity of water 

Friction coefficient 

Heat transfer coefficient 

Thermal conductivity of stainless steel 
Electrical power 


Velocity as measured by pitot tubes 


Average velocity 


Coefficients 

Effective inside diameter of test section 
Actual inside diameter of test section 
Outside diameter of test section 
Accelération of gravity 

Water mass flow rate 

Water pressure 

Heat transfer rate 

Bulk temperature of water 


Bulk temperature of water initially (i.e. 
at inlet to test section) 


Temperature of inside wall of test section 
Temperature of outside wall of test section 
Test section wall thickness 


Distance downstream from station No. l 
(see Fig. 3) 


éq in 


Btu/lb-°F 


Btu/ft-hr-CF 
watts 
ft/sec 


ft/sec 


in 

in 

in 
ft/sec“ 
lbs /sec 

psia 

Btu/sq in-sec 


ep 





Ap Pressure drop between stations noted in psi 


subscript 
At. Temperature drop through test section wall °F 
p Density of water lbs /f£t7 


subscripts <At stations No. 1, 2, 3 
lL, 2. 3 


i? 
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APPENDIX B 


Dimensions 


Test Section with Pitot Tubes Installed: 
Outside Diameter (d 3 
Wall Thickness (w,) 
Inside Diameter (d;) 
Length of Heating Section 


Heating Area (A), based on Inside 
Diameter 


Diameter of Pitot Tubes 


Test Section without Pitot Tubes: 
Outside Diameter (d o! 
Wall Thickness (w,) 

Inside-Diameter (d;) 

Lengtit.of Heating Section 


Heating Area (A), based on Inside 
Diameter 


0.3786 in 
0.0103 in 
0.3580 in 
17.75 in 


19.95 sq in 
0. 058 in 


0.3783 in 

0.0105 in 

0.3577 in 
17.75 in 


19.93 sqin 





sie 


APPENBIX C 


Sample Computations 


lor the sample cornputations, a particular run at a test sec 
tion pressure level of 265 psia, with a heat transfer (q/A) of 2.¢ 
Btu/sq in-sec, is here considered. The inlet bulk temperature (t,.) 
for this run was 307°F. The coefficients are computed for a point 
where the subcooling was 30°F. At the given pressure, 80°F sub- 


cooling corresponds to a bulk water temperature (t,) of 326°°. 


Friction Coefficient (C..): 
The pressure variation along the test section is assumed 


to be of the form, 
2 
psatbx +ex. (3) 


Using the experimental valuea at x = 0, 7 and 14 inches with the aid 
of Figs. 25 airl 26, the coefficients are obtained by means of the 


following three equations: 


265.90 


a 
264.35 =a+7b + 49c 
263.53 =a + 14b + 196c. 


When solved simultaneously, one obtains: 


zs 
iF 


-9. 0807, 


and 


-0.001735. 


ev) 
tt 
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The gradient, 
dp . 
<< =b + 2cx, (9) 


is found to be: 
SP = - (0.0807 + 0.00347x). 


The value of x (see Fig. 3) is determined by simple propor- 
tioning from a consideration of the actual bulk terznperature rise 
in the test section (see Fig. 8) and of the length of the test section 
(sce Appendix B). In this way, x is determined to be equal to 4.9 


inches. Then, 


se = 0.0977 lbs/in’, 


At the test section entrance the average velocity (¥V) of the 


ae 


water is given by: 


>. 
2 4m x 144 ; (6) 
px d' 
where d' is in this case the actual inside diameter (d.) of the test 
section. 


ds = d' 2 6.3580 in, 


i 


m = 0.73 lbs/sec 


and 
p = 57.0 ibs /it?, 
v P.. x0.73 x 144 
57. 0x3. 140.3585 


but V must be evaluated at x = 4.9 in. From Fig. 16 it can be seen 


= 138.3 {t/sec, 


that the velocity at station No. 2 must be corrected by the ratio of 
$3 . There is no observable increase in the measured velocity 





23%. 
at station No. 1, se no correction is made to the average velocity 
there. By linear interpolation to the value of x between the two sta- 
tions, 
V = 18.6 ft/sec. 


The friction coefficient is given by: 


1/2 
ad m 17239 : = 
Ce ae SITS TT (7) 
¥ p « 
Therefore, 
0.0977 = 32.2 x 0.73! /* x 1728 
Ce a oe = §.004] 


18.6 x 56.5 * 3.14 
Heat Transfer Coefficient (C,,): 
At. (actual) = 51 °x (see Fig. 8), 


A 219.95 in (see Appendix 5) 
C = 1.04 Btu/lb-°F. 
The effective heat transfer rate (q/A) is given by: 
’ mG 
= St. = - (12) 


= 1.94 Btu/sq in-sec. 


~Q 8 pP}a 


_ 51% 0.73 % 1.04 
a | ae 


The temperature of the outside wall (t_,) is obtained frorm 


Figs. 11 and 12 with linear interpolation between stations. 


t. 
wo 


The temperature drop (At) through the test section wali is given 


by (6): 
bee Bak (Ect) 200m es (3) 
~ 2 eae | 


WwW Oo 





eee. 
where q/A = 2.0 Btu/sq in-sec, since Eq. (3) was derived on the 
assumption of no heat loss to the surroundings. 


a, = 0.3786 in, 


and 
w, = 0.0103 in (see Appendix B). 
Therefore 
x oat 9.0103 [1 . 0.0103 3600 x 144 _ 33°R 
w 13.2 3 | Ox0. 12 al 
- t) = 475° ~ (33° + 326°) = L16°F. 
The heat transfer coefficient (C.) is given by: 
144 
CF ot (11) 
CH, 5-t) ,-t,) 
Therefore, 
144x 1.94 
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(Btu/sq inssec) 


Heat Transfer Rate 


Convection Ifucleate 


Boiling 


Wall Temperature (°F) 


Fige 1 Typical Heat Transfer Curve (1) 
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Fig. 3 Photograph of Pitot Tube Installed in 
Test Section 
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————> Test Cycle Flow Path 
——-—--—>  £«Heat Cycle Flow Path 


Throttle Valve 

Vent 

Receiver Tank 

Test Section 
Compressed Nitrogen 


Pressure Regulator 


Valve 
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Circulating Pump 
Pressure Gage 
Reservoir Tank 
iImnaersion Heater 
Venturi Flowmeter 
Thermoneter 


Fig. ) Test Circuit Diagran 
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Fig. 5 Photograph of Glass Sealed Thermocouple 
with Packing Gland 
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Fig. 6 Photograph of Test Section Installation 
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